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ABSTRACT 



o 

Context. Knowledge of radial velocities, elemental abundances, and accretion properties of members of star-forming 
regions is important for our understanding of stellar and planetary formation. While infrared observations reveal the 
evolutionary status of the disk, optical spectroscopy is fundamental to acquire information on the properties of the 
_ central star and on the accretion characteristics. 

Aims. Existing 2MASS archive data and the Spitzer c2d survey of the Chamaeleon II dark cloud have provided disk 
properties of a large number of young stars. We complement these data with optical spectroscopy with the aim of 
providing physical stellar parameters and accretion properties. 

Methods. We use FLAMES/UVES and FLAMES/GIRAFFE spectroscopic observations of 40 members of the 
Chamaeleon II star-forming region to measure radial velocities through cross-correlation technique, lithium abundances 
by means of curveso of growth, and for a suitable star elemental abundances of Fe, Al, Si, Ca, Ti, and Ni using the code 
MOOG. From the equivalent widths of the Ha, H/3, and the He I A5876, A6678, A7065 A emission lines, we estimate 
C*| . the mass accretion rates, M acc , for all the objects. 

Results. We derive a radial velocity distribution for the Chamaeleon II stars, which is peaked at (K a d) = 11.4±2.0 km s" 1 . 
We find dependencies of M acc oc M*' 3 and of M acc oc Age~ 082 in the ~ 0.1 — I.OMq mass regime, as well as a mean 
mass accretion rate for Chamaeleon II of M acc ~ 7lg 6 x 1O~ 1O M0 yr" 1 . We also establish a relationship between the 
He i A7065 A line emission and the accretion luminosity. 

Conclusions. The radial velocity distributions of stars and gas in Chamaeleon II are consistent. The spread in M acc at 
a given stellar mass is about one order of magnitude and can not be ascribed entirely to short timescale variability. 
Analyzing the relation between M acc and the colors in Spitzer c2d and 2MASS bands, we find indications that the inner 
disk changes from optically thick to optically thin at Macc ~ 10~ 10 Af© yr -1 . Finally, the disk fraction is consistent with 
the age of Chamaeleon II. 

Key words. Accretion - Stars: pre-main sequence/low-mass/abundances - Open clusters and associations: individual: 
(f) ■ Chamaeleon II - Techniques: spectroscopic 
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CN| ■ 1201 lh . The rate at which the central star accretes from disk 
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The study of accretion properties of members of star- mate rial has been found to approximately scale with the 
forming regions (SFRs) is important for our understand- squ a re of the stellar mass and to decrease with age (see, 
ing of stellar and planetary formation. While infrared ob- e g ; iHerczeg k Hillenbrand! I2008L and references therein), 
servations provide information on the structure of the cir- In addition, the accretion properties are also important to 
?J ; cumstellar disk, the accretion properties can be retrieved understand the planet-metallicity relation. In fact, the effi- 
CS . from photometry and spectroscopy using primary diagnos- ciency of dispersal of circumstellar disks is predicted to de- 
fies, suchas^hej^excessemission (e.g., |Gullbring et al.| perid on stellar me tallicity in the sense that the formation 
[1223 IRigliaco et al.| |2011a|), the Paschen/Balmer con- of planetesimals around star s is faster at higher metallic- 
tinuum and Balmer lump (e.g., |Gullbring et_ajj ll998t ity dErcolano fc Clarkdl2lHflh . Simultaneous measurements 
IHerczeg k Hillenbrand 2008; Rigliaco et al.||2011b|), or sec- of accretion rates and elemental abundances in SFRs and 



ondary tracers, like hydrogen recombination lines (Ha, y0U ng clusters are therefore crucial to shed light on the role 

H/3, H 7 , H9, Pa/3, Pa 7 , BrA and the Hel, Call. Nal of metallicity in disk dispersal and planetary formation, 
lines (e.g., iMuzerolle etaH Il998alfbfc iNatta et al.l 120061: 

The Chamaeleon II ( hereafter Cha II) da rk cloud, at a 

~~Z j ~Z "~ , , jr n . distance of 178 ± 18 pc (IWhittet et alJll997fl . is one of the 

bend offprint requests to: K. Biazzo . , . . 1 „. — — : — , „>. 

* ti j i-iT /^TnAnnn , TT^^T^o^ u j.- i three mam clouds ol the Chamaeleon complex (a ~ 12 , 

Based on FLAMES (GIRAFFE+UVES) observations col- 2 1 

lected at the Very Large Telescope (VLT; Paranal, Chile). 5 ~ - 78 ). It extends over ~ 2 deg in the sky fsee lLuhmanl 

Program 076.C-0385(A). 2008 for a recent review). The population of Cha II consists 

Correspondence to: katia.biazzo@oacn.inaf.it of some 20 classical T Tauri stars (CTTSs), ~ten weak- 
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lined T Tauri stars (WTTSs), an inter mediate-mass Herbig 
Ae st ar (IRAS 12496—7650; see, e .g.. iGarcia Lopez et al.l 
1201 If) , a few Herbig-Haro objects (|Alcala et al.l 120081 and 
references therein), ~three sub-stellar object s, and refer- 
ences therein) and five very low-mass stars (jSpezzi et al.l 
120081 and references therein) . Cha II is one of the five SFRs 
included in the Spitzer Space Telescope Legacy Program 
"From Molecular Cores to Planet-f o rming Disks" (c2d \ 
lEvans et al.l 120031: lYoung et al.l 120051: iPorras et all [2007). 
Through extensive work based on c2d IRAC and MIPS 
Spitzer fluxes and complementary data, a reliable cen- 
sus of the population in Cha II (down to 0.03A/q) was 
achieved by lAlcala et al.l (J2008) . They concluded that the 
cloud is dominated by objects with active accretion, with 
the Class II sources representing ~ 60%. The same sam- 
ple was investigated spectrosco pically in th e optic al us- 
ing GIRAFFE/UVES@VLT by ISpezzi et~aH (|200l . who 
derived stellar parameters and estimated a mean age of 
4 ± 2 Myr for Cha II. However, studies of radial veloci- 
ties, elemental abundances, and accretion properties of the 
cloud members were not addressed in these works. Recently, 
lAlcala et all (|2011al) analyzed the star IRAS 12556-7731, 
concluding that it is indeed a background lithium-rich M- 
giant star unrelated to Cha II. 

As a continua t ion of the studies by lAlcala et al.l (|2008h 
and ISpezzi et alj (|2008l ). we derive here radial velocities, 
elemental abundances, and accretion properties for 40 pre- 
main sequence (PMS) stars in Cha II. The outline of the 
paper is as follows. In Sect. 03 we describe the spectro- 
scopic observations and data reduction. In Sect. [21 we re- 
port determinations of radial velocities, elemental abun- 
dances, and accretion properties. The main results are dis- 
cussed in Sect. 01 while our conclusions are presented in 
Sect. [5] 



2. Spectroscopic observations and data reduction 

The observations were conducted in February-March 2006 
and February 2007 using FLAMES (GIRAFFE+UVES) at 
the VLT. A complete jou rnal of the ob s ervati ons and instru- 
mental setup is given in ISpezzi et al.l (|2008f ) . The relevant 
information for this paper is summarized in Table [T] 

We observed 32 objects with GIRAFFE, 11 with UVES, 
and two with both spectrographs (see Table [lj. Despite the 
~ 25' FLAMES field of view (FoV), it was not possible to 
assign a large number of fibers in each configuration to the 
PMS objects/candidates because of the large spatial scat- 
ter of the targets. The remaining fibers were allocated to 
young candidates and field stars_l (nu mbers in parenthe- 
ses in columns 4 and 5 of Table [TJ see ISpezzi et al.ll2008l 
for more details). Thirty- two objects were observed several 
(2-4) tim es within two days ( see Table 03). 

While ISpezzi et all (H0O8) used a single spectrum per 
object to derive the spectral type and to confirm the 
presence of Lil 6708 A absorption, we use here the com- 
plete set of spectra to investigate accretion and short 
timescale variability. To this aim, we reprocessed the 
FLAMES/GIRAFFE and FLAMES/UVES observations. 
The GIRAFFE data were reduced using the GIRAFFE 
Base-Line Data R eduction Software 1.13.1 (girBLDRS; 
IBlecha et aLll2000h : bias and flat subtraction, correction 



1 Hereafter, we refer to this sample as "field stars" . 



Table 1. Summary of the observations. 



Instrument 


Range 


Resolution # stars 


# spectra 




(A) 


(A/AA) 




GIRAFFE 


6438-7164 


8 600 32(+27) 


69(+45) 


UVES 


4764-6820 


47000 11 


25 



for the fiber transmission coefficient, wavelength calibra- 
tion, and science frame extraction were performed. Then, 
a sky correction was applied to each stellar spectrum using 
the task SARITH in the IRAP0 echelle package and by 
subtracting the average of several sky spectra obtained si- 
multaneously. The reduction of the UVES spectra was per- 
forme d using the pipeline developed by iModigliani et al.l 
(|2004l) . which includes the following steps: subtraction of a 
master bias, echelle order definition, extraction of thorium- 
argon spectra, normalization of a master flat-field, frame 
extraction, wavelength calibration, and correction of the 
science frame for the normalized master flat-field. Sky sub- 
traction was also performed with the IRAF task SARITH 
using the fibers allocated to the sky. 

3. Data analysis and results 

3.1. Radial velocity distribution, membership, and binarity 

We determined radial velocities (RVs) of each object, choos- 
ing Hn 23 and RX J1303.1-7706 as UVES and GIRAFFE 
templates, respectively. These slowly rotating stars show 
no strong accretion signatures (see Table 03) . We measured 
the RV of each template (highest S/N) spectrum using the 
IRAF task rvidlines inside the RV package, which consid- 
ers a line list. We used 50 and 10 lines for the UVES and 
GIRAFFE spectra, respectively, obtaining V ra( j = 12.5 ±0.4 
km s- 1 for RX J1303. 1-7706, and V rad = 15.2±0.3 km s" 1 
for Hn 23. The heliocentric RV of all targets was determined 
through the task FXCOR of the IRAF package RV, which 
cross-correlates the target and template spectra, excluding 
regions affected by broad lines or prominent telluric fea- 
tures. The centroids of the cross-correlation function (CCF) 
peaks were determined by adopting Gaussian fits, and the 
RV errors were computed using a procedure that consid- 
ers the fitted peak he ight and the antisymmetric noise (see 
iTonrv fc Davislll979h . When more spectra were acquired, 
we computed the average RV for each object. 

In order to estimate the binary fraction, we considered 
as singles the stars with only one CCF peak and with night- 
to-night RV variations within 3cr. In the end, excluding 
seven objects for which we could not measure the RV, we 
find six spectroscopic binaries, which means a binary frac- 
tion of 18%. In the last column of Table 03 we list the most 
probable spectroscopic systems. Five stars are newly identi- 
fied as spectroscopic binaries, while RXJ1301.0— 7654a was 
already recognized as a do uble-lined (SB2) spectroscopic 
system (|Covino et al.lll997h . 

In Fig. [Q we show the Cha II distribution of the RV 
measurements in the local standard of rest (LSR) obtained 

2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of the 
Universities for Research in Astronomy, inc. (AURA) under co- 
operative agreement with the National Science Foundation. 



2 



K. Biazzo et al.: Radial velocity, abundances, and accretion of Chamaeleon II 



from both the UVES and GIRAFFE s pectra, along w i th the 
RV distribution of the gas derived bv lMizuno et al.l (|l999) 
from the C 18 ( J = 1 — 0) transition in 11 dense molec- 
ular cores in Cha II. Excluding the spectroscopic binaries, 
the RV distribution of the Cha II members has a mean of 
(V ra d) = 11 A ± 2.0 km s _1 , which translates to a value of 
(Vlsr) = 4.9 ± 2.0 km s" 1 . The Gaussian fit of the distri- 
bution yields a mean value peaked at (Vlsr) = 3.9 ±1.6 
km s~\ i.e. (V-ad) = 10.0 ± 2.1 km s _1 , which is m 
fairly good agreement with the average velocity of the gas 
((Vlsr) = 3.0 ± 0.7 km s"\ i.e. (V rad ) = 9.6 ± 0.7 km s" 1 
jMizuno et al.lll999t) . 



This work 

Mizuno et al, (1999) 




20 

(km/ s) 



30 



40 



Fig. 1. Average RV distribution in the LSR (solid thick 
line) of the Cha II PMS s tars. T he distribution of the gas 
derived by iMizuno et al.l (|1999t) is overlaid (dashed line) . 
The Gaussian fit to the PMS RV distribution is shown (thin 
line). The shaded histogram represents spectroscopic binary 
stars, while the hatched one marks the UVES observations. 
In the case of Hn 24 and Sz 54, where both UVES and 
GIRAFFE RVs were measured, we considered the UVES 
observations. 



(RXJ1301.0— 7654a and Sz 54), while the rest can be con- 
sidered as field stars. For the SB2 stars, the average of the 
RVs of their components falls inside the (Vr a d) ± 3<r distri- 
bution. The relatively narrow RV distribution of the strong- 
lithium sample confirms that these stars are all members 
of the same association. 



0.6 



0.4 - 



0.7 




1.0 - X 



Fig. 2. EWu versus RV for stars in the Cha II FoV (see 
Table [JJ). Open circles refer to most probable single stars, 
filled circles are multiple components, and asterisks refer to 
field stars. We excluded stars whose binarity was detected 
from RV variation at different phases (namely, Sz51 and 
Hn 24; see Table[S|). In the case of Sz54, observed with both 
FLAMES configurations, we considered only the UVES RV 
values. Vertical lines represent the (V-ad) ±3cr values, where 
cr=2.0 km s _1 (see Section u3~T|) . 



3.2. Lithium equivalent width and radial velocity 

Lithium equivalent widths {EWu) were measured by di- 
rect integration or by Gaussian fit using the IRAF task 
SPLOT. Errors in EWu were estimated in the following 
way: i) when only one spectrum was available, the stan- 
dard deviation of three EWu measurements was adopted; 
ii) when more than one spectrum was gathered, the stan- 
dard deviation of the measurements on the different spectra 
was adopted. Typical errors in EWu are 

of 0.001-0.087 A 

(higher values for GIRAFFE data); for the stars C62 and 
C66, (Tew u ~0.15A. Our EWu measurements are consis- 
tent with the values of ISpezzi etaH (|2008l ) within 0.02 A. 

Figure [2] shows the EWu versus RV for the 37 Cha II 
members listed in Table [5] (circles) and for the field stars 
(asterisks). The difference in RV distribution between the 
no-lithium or weak-lithium stars, and the strong-lithium 
stars is noticeable. The strong Li stars are confined to a 
narrow range of velocities (i.e., inside ±3<r from the peak of 
the RV distribution). The strong-lithium sample contains 
single stars and also the two most probable SB1 systems 
(C66 and IRAS F13052-7653N) and the two SB2 systems 



3.3. Elemental abundances 
3.3.1. Abundance measurements 

The FLAMES/UVES wide spectral coverage allows us to 
select several tens of Fei±Feii lines and spectral features 
of other elements to m easure abundances f rom line EWs. 
To this aim, as done in lBiazzo et al.l (|2011bj ). we discarded 
stars with T c g S 4000 K (because of significant formation 
of molecules in the atmosphere), fast rotators (to avoid ro- 
tational blending) , and strong accretors (for which accurate 
abundance analysis is hampered) . In the end, only one star 
(Hn 23) fulfills the required criteria. Effective temperature 
and surface gravity (log g) from the literature (Table [7]) 
were used as initial values, and initial microturbulcnce (£) 
was set to 1.5 kms -1 . Final values of the atmospheric pa- 
rameters are listed in Tabic [2j together with abundance 
determinations, abundance internal errors, and number of 
lines considered (in parenthesis). The first source of inter- 
nal error in abundance is due to uncertainties in line EWs, 
while another contribution comes from the uncertainties in 
stellar parameters. Systematic (external) errors, introduced 
by the code and/or model atmospher e, are negligible in 
comparison with the internal ones (see[B lazzo 
for details on the treatment of errors). 
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Table 2. Spectroscopic parameters and elemental abun- 
dances of Hn 23. 



Spectroscopic parameters 

T efi (K) 4500 ± 100 
logg (dex) 4.0 ±0.2 

£ (kms' 1 ) 1.9 ±0.2 

Elemental abundances 



[Fel/H] 

[Fen/H] 

[Al/Fe] 

[Si/Fe] 

[Ca/Fe] 

[Tii/Fe] 

[Tin/Fe] 

[Ni/Fe] 



-0.12±0. 
-0.13±0. 
-0.02±0. 
±0.02±0. 
-0.07±0. 
-0.20±0. 
±0.04±0. 
±0.01±0. 



14±0 
21±0 
14±0 
15±0 
15±0 
18±0 
16±0 
16±0 



06(41) 

06(3) 

05(2) 

07(1) 

10(3) 

09(5) 

08(1) 

04(7) 



The [Fe/H] value of Hn 23 is slightly below the so- 
lar value and in agreement with the mean abund a nce o f 
<[Fe/H]>= -0.11 ± 0.11 found by ISantos et all (l2008h 
for the Chamaeleon Complex (see Fig. |3J). This supports 
the suggestion that SFRs in the solar neighborhood are 
slightly more metal- poor than nearby young open clusters 
(Bia zzo et all 1201 lal ). This issue certainly deserves further 
study, using more stars in the region and homogeneous sam- 
ples of PMS stars in as many SFRs as possible. Moreover, 
all other [X/Fe] abundances are close to the solar ones, 
with silicon and nickel close to the cluster mean value of 
<[ Si/Fel>= 0.03 ±0.01 and <[Ni/Fe]>= -0.05±0.02 found 
bv ISantos et~aTI ([2008). Titanium seems to be affected by 
NLTE effects (see Tabled , as previously found by other au- 
thor s for stars with temperatures cooler than ~ 500 K (see, 
e.g.. lD'Orazi fc Randichll2009HBiazzo et al.ll2011al and ref- 
erences therein). However, detailed treatment of NLTE ef- 
fects is beyond the scope of this paper. 



0.4 - ■ This work 
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Fig. 3. Iron, silicon, and nickel abundances versus spectro- 
scopic temperature for Hn 23 and the stars analyzed by 
ISantos et ail (|2008l ). 



logn(Li) values). The average of logn(Li) is about 2.5 dex 
with a dispersion of 0.6 dex. The lowest lithium abundance 
values are presumably due to spectral veiling, which affects 
the line EW. 



3.3.2. Lithium abundance 

Mean lithium abundances were estimated fro m the average 
EWu listed in Table [5] and T c g values from ISpezzi et al.l 
(200 S), by using the LTE curves-of-growth reported by 
Pavlenko fc Maeazzul (|l996h for T eff > 3500 K, and by 
Palla et al.l (|2007t ) for T eS < 3500 K. The log g values were 
derived using the effective temperat ure, luminosi t y, an d 
mean mass reported for each star in ISpezzi et al.l (|2008j) . 
The main source of error in logrt(Li) co mes from the uncer - 
tainty in T cff , which is AT off - 100 K (jSpezzi et al.ll2008l) . 
Taking this value and a mean error of 0.020 A in EWu 
into account, we estimate a mean log n(Li) error ranging 
from -0.07-0.10 dex for cooler stars (T eff - 3200 K) down 
to ~0.05-0.09 dex for warmer stars (T cff ~ 4200 K), de- 
pending on the EWu value. Moreover, the logg value af- 
fects the lithium abundance, in the sense that the lower 
the surface gravity the higher the lithium abundance, and 
vice versa. In particular, the difference in logn(Li) may 
rise to —±0.05 dex when considering stars with mean val- 
ues of EW U = 0.500 A and T off = 4000 K and assuming 
A log g = ±0.5 dex. 

In Fig. 2] we show the mean lithium abundance as a 
function of the effective temperature (see Table \5\ for the 



3.4. Accretion diagnostics and mass accretion rates 

The spectral coverage of our data allows us to select sev- 
eral lines (namely, Ha A6563 A, H/3 A4861 A, Hei A5876 
A, He i A6678 A, and Hei A7065 A) that can be used 
to determine the accretion luminosity (L x cc ). These emis- 
sion lines are powered by processes related to accretio n 
from the circumstellar disk (jHerczee fc H illenbrand 2008). 
The use of these lines as secondary accretion diagnos- 
tics relies on empirical linear relationships between the 
observed line luminosity (L x ) and th e accretion luminos- 
ity (e.g., lHerczeg fc Hillenbr and 2008]) ■ These relationships 
have been established t hrough primary diagno stics, such 
as UV excess emission (jGullbring et al.l Il998l). We used 
the L x — L x cc empirical relations of iHerczeg fc Hillenbrand! 
(2008) to derive L x cc . The line luminosity was calculated 
as L A = ATrR'jF x : where the stellar radius, i?*, was taken 
from ISpezzi et al.l (|2008l) and the observed flux at the stel- 
lar radius, F x , was derived by multiplying the EW of each 
line (EW\) by the continuum flux at wavelengths adjacent 
to the line (-F C o„ t 1nuum)- The l atter was gathered from the 
Next Gen Model Atmospheres (jHauschildt et all 1 19991 ). as- 
suming the corresponding stellar temperature and gravity 
(see Table [7]). The mass accretion rate, M x cc , was then de- 
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Fig. 4. Upper panel: EWu versus effective temperature. 
The upper enve l ope fo r the Pleiades, as adapted by 
ISoderblom et al.l (|1993[ ). is overplotted as a dotted line. 
Lower panel: Lithium abundance versus effective temper- 
ature. The "lithium isochrones" by iD'Antona fc Mazzitellil 
(|1997|) in the 2-20 Myr range are overlaid with dashed lines. 
In both panels, empty symbols represent spectroscopic bi- 
naries, while mean error bars are overplotted on the lower 
right corner. 



rived from L^ cc using the following relationship (|Hartmannl 
1998): 



GM* 



(1) 



where the s tellar radius R+ and mass M+ for each star were 
taken from ISpezzi et al.l (12008T). and the inner -disk radius 
i?i n , when available, from lAkala etal] |2003)- When_ no 
R ln w as available, we assumed R ln = 5i?* (see iHartmannl 
119981) . which is a good approximation for most accretors, 
as pointed out by lAlcala et al.l ()2011bf) . Contributions to 
the error budget on M acc include uncertainties on stellar 
mass, stellar radius, inner-disk radi us, and L^ nn . Assum ing 
mean errors of ~ 0.15Mrp i n M* (|Spezzi et all 120081) . ~ 
O-IORq in TL (ISpezzi et alj I2008D. and - 0.2 AU in R in 
(jAlcala et al.ll2008h . 5 - 10% as relative error in EW X , 10% 
in F„_i7^ , and the uncerta inties in the relationships by 



iHerczee fc Hillenbrand! (j2008T ) . we estimate a typical error 
in log Macc °f ~ 0.5 dex. 

Apart from variability phenomena, which will be dis- 
cussed in Section 14.21 the mass accretion rates derived 
from the various diagnostics should be consistent with each 
other. In the following, we describe the results drawn from 
the hydrogen and helium emission lines. 



Figure [5] (upper panel) shows the mean accretion lu- 
minositjQ, derived from the Ha emission line, versus the 
mean accretion luminosities obtained from other hydrogen 
and helium emission lines. A fairly good agreement between 
these diagnostics is present. Comparing the mean mass ac- 
cretion rate from the Ha line with the mean M acc as ob- 
tained through the other diagnostics, the agreement is well 
reproduced (see lower panel in Fig. [5] and column 16 in 
Table [5]). This justifies the use of all the diagnostics to 
compute an average (L acc )EI (and, hence, also an average 
(Ma.cc)) for each star. A weighted average of (L acc ) derived 
from the Ha, H/3, Hei A5876 A, and Hei A6678 A emission 
lines allows us to analyze the relationship between (L acc ) 
and the luminosity in the He I A7065 A line, which is shown 
in Fig. [5J A linear fit to the relationship gives 



(logL 



HcA7065\ 
acc / 



(4.01±0.35)+(1.14±0.06)(logL 



HoA7065 



>-(2) 



The good correlation justifies the use of the He I A7065 A 
line as the reliable diagnostic of L acc - 

As already pointed out, the different line diagnostics 
yield consistent mass accretion rates (see Fig. [5] and Table 
E]). Considering, for instance, the Ha line, which is observed 
in all targets, the mean difference in log M^ c " as compared 
to log-M^cc* is of 0.3 ± 0.3 M Q yr _1 (with a maximum of 
0.8 dex observed for Sz 54), while it is -0.2 ± 0.4 dex 
with respect to logM^ A5876 (with a maximum of —0.2 dex 
observed for Sz 56), and 0.1 ± 0.5 dex with respect to 
logM^f 6678 (with a maximum of 1.1 dex for Sz 50). 
The mass accretion rate for the sample is in the range 
10~ n ~2 x 10~ 8 MQyr" 1 , which is typical of Class I I 
low-mass YSOs (see Fig. 2 in ISicilia-Aguilar et al.ll2010h . 
Excluding multiple systems, stars with EWna < 10 A, and 
the early-type star DK Cha, we find an average mass ac- 
cretion rate for Cha II of (M aGC ) ~7x 10~ 10 M yr _1 . 

3.4.1. Comparison with previous studies 

Fourteen of o ur tar gets were also analyzed by 
lAntoniucci et all (|2011h as part of the POISSON 
(Protostellar Objects IR-optical Spectral Survey On 
NTT) project aimed at deriving the mass accretion 
rates of young stars in star-forming regi ons through low- 
resolu tion optical/near-IR spectroscopy. lAntoniucci et al.l 
(|201lD used the Br7 line as an accretion tracer. They 
argue that this is the best diagnostic in their spectra when 
compared with other tracers (i.e., Pa7, Can, Ha, and 
[Oi]). Comparing their logM acc with our mean logM acc , 
their values tend to be systematically higher than ours. 
In particular, their values strongly diverge from ours at 
logMacc < —8, with a mean difference of 1.0 ± 0.8 dex. 
An analogous trend was found by the same authors when 
comparing ^ ac i c 9 ^ 0l ^ HQ,Can with £ a< 5<J ( see their Fig. 4), 
which is more evident in their i. aC c — ^acc diagram. 
They ascribe this behavior to enhanced chromospheric 
emission or absorption from outflowing material, direct 
photoionization at higher luminosities, or flux losses due 
to winds. 

When comparing the POISSON's EWn a with the val- 
ues derived by us (left-hand panel of Fig. [7]), the difference 



3 This is the average of different observations for a given star. 

4 This is the average of different diagnostics. 
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Fig. 5. Average L^ c " (top panel) and M^™ (bottom panel) 
derived from the Ha line as a function of (£ acc ) and (M acc ) 
obtained from the H/3, He I A5876 A, and He I A6678 A lines. 
Mean error bars are overplotted on the lower right corner 
of each panel. 



is AEWn a = 10 ± 23 A. Considering that the observa- 
tions were performed at different epochs (our run was in 
2006, while their run was in 2009) and that their EW^a 
were measured on low-resolution spectra (R ~ 700), the 
agreement is fairly good. In fact, a good correlation is 
found between our estimated Ha line luminosities and the 
POISSON's values, with a difference of 0.04 ± 0.13L© (see 
Fig. middle panel). Thus, the differences in mass accre- 
tion rate arise when deriving L acc . A similar behavior as 
the one shown in Fig.[7]for the mass accretion rate is found 
when our L acc values are compare d with those of PO I SSON , 
in agreement with the results bv lAntoniucci et al.l (|201l|) . 
This means that the differences between POISSON and our 
determinations are mainly due to the tracers used to derive 
-L acc . As shown in Fig. [5J very good correlations of L acc 



as derived from the Hei and Ha lines are found. Should 
there be an important contribution to the line diagnostics 
by winds and/or chromospheric activity, our estimates of 
L aC c would be in excess with respect to those derived from 
the Br7 line. Instead, the opposite is observed. We thus 
exclude the possibility that the differences are due to the 
influence of winds and chromospheric activity. A possible 
explanation for the difference between our average L acc val- 
ues and the POISSON's may be that the various ac- 
cretion diagnostics originate in different regions at different 
physical conditions. In order to investigate this issue, mea- 
surements of primary accretion tracers, and Br7 measure- 
ments are needed. Such analysis cannot be conducted with 
the data available here, but will be addressed in our future 
studie s exploiting X-Shooter@VLT data (cf. lAlcala et al.l 
[2011hh . 







A 



o 



He I A7065 



/ 

/ 



4 
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log (L He706 VL Q ) 
Fig. 6. Average accretion luminosity versus line luminosity 
for the He I A7065A line. The linear fit given in the text 
(Equation [2]) is represented by the dashed line. Mean error 
bars are overplotted on the lower right corner. 



3.4.2. Other lines in the optical 

A number of optical lines are seen in emission in the spec- 
tra of several stars in our sample. These are the forbid- 
den lines of [Oi] A6300.3 A, and A6363.8A, [Sn] A6715.8A, 
and A6729.8A, [Nn] A6548.4A, and A6583.4A, as well as 
Fell multiplets. These lines trace mainly stellar /disk winds , 
jets, disk surfaces, an d outflow activity (jCabrit et al.l fl990; 
lHartigan et al.lfl995fl . and their detailed treatment is be- 
yond the scope of this paper. Here, we only note that 
the [Oi] A6300.3A line is detected in 6/11 objects, the 
[Oi] A6363.8A line in 1/11 objects, the [Sn] 6715.8A and 
6729.8 A emission is observed in 7/40 and 9/40 sources, 
respectively, and the [Nil] 6548.4 A and 6583.4 A lines 
are present in 5/40 and 5 /40 stars (see Table [5]) . In the 
case of Sz51, the star showing the strongest Ha in our 
FLAMES/UVES sample, there is evidence of other emis- 
sion lines (such as the Mgi triplet at A5167.3, 5172.7, 
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5183.6 A, the multiplets 42 and 49 of Fen, etc.) indicat- 
ing mass loss. 



4. Discussion 

4.1. Accretion versus stellar age and mass 

Figure \8\ shows the mean mass accretion rate versus stel- 
lar age (see Table [7]) for all the targets except the early- 
type star DK Cha, which ha s a massive disk and high M acc 
(jSicilia-Aguilar et all [2010). A slightly decreasing trend 
with age may be present, though over a narrow age in- 
terval (~ 0.4 — 13.4 Myr) and with a large scatter in M acc . 
This would be consistent with the evolution of a viscous 
disk (see, e.g.. lHartmann et al.lll998l ISicilia-Aguilar et al.l 
|2010L and references therein), although the mean M acc 
is slightly lower than that expected from the model at 
the Cha II age. In order to quantify the degree of anti- 
correlation between log M acc and log Age, we calculated the 
Spearman's rank correlat ion coefficient usi ng the IDL pro- 
cedure R_CORRELATE ([Press et al.lll996ft . We find a cor- 
relation coefficient p = —0.38, with a probability of obtain- 
ing such p from randomly distributed data of p = 0.08. 
This seems to confirm a moderate anti-correlation between 
Macc and age, with an average mass accretion rate of 




-0.5 



0.0 0.5 

< log Age> (Myr) 



7lf x 10 _10 M Q yr- 1 at a mean age of Myr. The lin 
ear relation we obtain for the likely single stars with mean 
EW}i a higher than 10 A is 

(logM acc (Moyr^ 1 )) = -8.84±0.22-0.82±0.40 logAge (Myr 

The slope is higher than, yet consistent with i n the er- 
rors with, that obtained by lHartmann et al.l (|1998|) in 
Cha I (i.e., logM acc (Moyr" 1 ) = -8.00 ± 0.10 - 1.40 ± 
0.29 log Age (Myr)). We warn, however, about the large un- 
certainties in absolute ages derived from theor etical models 
for s tars younger than ~ 10 Myr (see, e.g., ISpezzi et al.l 
2008). We note also that a strong constraint on the appar- 
ent trend in logM acc versus age is set by only one object 
in the sample (C41, age ~ 13 Myr). The trend disappears 
if this object is not considered. 

Figure [S] shows the mean mass accretion rate versus 
mean stellar mass (see Table [JJ . Taking into account sin- 



Fig. 8. Mean mass accretion rate versus age. Open sym- 
bols represent the targets with mean EWn a < 10 A, while 
binaries are evidenced as upper limits. The big square rep- 
resents the mean position of Cha II, taking into account 
the single stars with mean EWn a > 10 A (vertical and hor- 
izontal error bars correspond to the standard deviations 
from the mean M acc and age, respectively). The dashed 
lines mark minimum and maximum limits of Eq. [3j while 
dotted lines represent th e analogous limits of th e relation- 
ship derived for Cha I by lHartmann et al.l ([1998) ■ The col- 
lection of viscous disk evolutionary models for solar-type 
stars with initial disk masses of 0.1 — 0.2M©, constant vis- 
cosity a = 10~ 2 , and visco sity exponent 7 = 1 reported by 
ISicilia-Aguilar et all (|2010h is also displayed (filled region) . 



gle stars with EWn a > 10 A, we find the following linear 
relationship: 

(logM acc (Moyr- 1 )) = -8.64±0.21+1.30±0.41 logM* (M ) ,(4) 

where th e slope is slightly lower than that found for Cha I 
(<~ 2; see lAntoniucci et alJl20lil and references therein). A 
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maximum limits of Eq. |4j while the dott ed line represents 
the m ean relationship derived for Cha I bv lAntoniucci et"af1 
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corner. 



correlation of M acc with stellar mass is evident. Spearman's 
rank correlation coefficient is p — 0.51 with a probability 
p = 0.02. The exponent of the M acc — M* power law is con- 
sistent with the range of values ~ 1.0 — 2.8 found for low- 
mass stars in other SFRs (see, e.g., Herczeg k, Hillenbrand! 
120081 iFang et al.l 120091: lAntoniucci et all 1201 ll and refer- 
ences therein). 



4.2. Short timescale variability 

Young PMS stars are known to be variable, due to the 
combination of different processes dHerbst et al.lll994h oc- 
curring on different timescales: on short timescale (days), 
variability can be induced by rotation of cool or hot spots 
(type I variability), and on long timescales (months-years), 
accretion rate changes (type II variability) or obscuration 
by circumstellar dus t might occur (type III variability, e.g., 
Schi sano et al.ll2009h . The time span of our observations is 
~ 48 hours. Emission line variability on timescales shorter 
than two days is observed for several objects and for all the 
ana lyzed lines (see Fig. IB.1[) . In particular, Figs. IA.l1 IA.21 
and lA.3l show the Ha and H/3 profiles of the stars displaying 
the strongest line variations, while Fig. IB.ll shows the mass 
accretion rates derived from the Ha, H/3, He I A5876 A, and 
He i A6678 A lines versus stellar mass. The vertical bars in 
these plots represent the range of M acc as due to the two- 
day variability of the corresponding diagnostics. 

In general, on the timescale of only two days, the Ha 
equivalent width of some stars changes up to a factor of 
2 — 3 and, assuming that this is due to accretion, variation 
in log M acc would be 0.2 — 0.6 dex, i.e., a factor of 1.6 — 4.0. 
We conclude that this logM acc variability, even if large, 
cannot explain the logM acc spread at a given mass. This 



means that other stellar properties besides mass must also 
affect the variations. 



4.3. Accretion and Metallicity 

Investigating the dependency of the mass accretion rate 
upon iron abundance in SFRs is important for two rea- 
sons. First, while the correlation between stellar metallic- 
ity and presence of giant pl anets around solar-ty pe stars 
is well established (see, e.g., iJohnson et all [20101 and ref- 
erences therein), the metallicity-planet connection in the 
early stages of planetary formation is still a matter of de- 
bate. The evolution of M acc is affected by possible planetary 
formation in the disk, and hence it might provide important 
clues on the planet-metallicity correlation. Second, the effi- 
ciency of the dispersal of circumstellar (or protoplanetary) 
disks and hence the dispersal timescale are predicted to de- 
pend on metallicity in the sense that pla netary formation is 
faster in disks wit h high er metallicity ([Ercolano fc Clarke] 
l2010h . lYasui et all <|2010h find that the disk fraction in sig- 
nificantly low-metallicity clusters ([0/H]~ —0.7) declines 
rapidly in < 1 Myr, which is much faster than the value of 
^5 — 7 Myr observed in solar-metallicity clusters. Since the 
shorter disk lifetime reduces the time available for plane- 
tary formation, they suggest that this could be one of the 
reasons for the strong pla net-metallicity cor rel ation - 
Recent s tudies by iDe Marchi et all (|201ll ) and 
iSpezzi et all (|2012[) in the Large and Small Magellanic 
Clouds show that metal-poor stars accrete at higher rates 
compared with solar-metallicity stars in galactic SFRs. 
Summarizing the mean M acc of low-mass (0.1 — I.OMq) 
Class II stars members of ~ 3 — 4 Myr old nearby SFRs 
for which iron abundance has been recently measured (see 
Table 13]), it is only possible to point out that for [Fe/H]~ 
the mass accretion rate is ~ 10 -10 Mq yr -1 . 



4.4. Fraction of accretors versus age 

Excluding the six spectroscopic binaries, 27 of the studied 
stars result in having mean EWb_o higher than 10 A, which 
would imply a percentage of accretors of about 26/34=79% 
(34 being the total number of single stars in the sample). 
However, the majority of PMS stars in Cha II have a spec- 
tral type later than K7, with most of t hem later than 
M3. T herefore, according to the criterion by I White &: Basrl 
«2003h . a more adequate dividing line between most prob- 
able accretors and non-accretors in Cha II is EWn a =20 A. 
Using this criterion, 19 stars can be classified as true accre- 
tors, leading to a percentage of ~ 55 ± 5%. This fraction 
of accretors is consistent with the average age of the cloud 
members. In fact, following the m ass accretion rate e volu- 
tion with time shown in Fig. 3 of iFedele et al.l (|2010l ). the 
fraction of stars with ongoing mass accretion decreases fast 
with age, going from ~ 60% at 1.5 — 2.0 Myr down to ~ 2% 
at 10 Myr. 



4.5. Color-M aC c diagrams 

Near-in frared colors ca n be u sed to probe the inner disk 
region. lHartigan et al.l (jl995l) . studying a sample of 42 
T Tauri stars and using the K — L color excesses, pointed 
out that the disk dispersion is mainly due to the for- 
mation of micron-sized dust particles, which combine to 
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Table 3. Iron abundances and mass accretion rates in nearby SFRs. 



Star- Forming Region P e /H] Reference 

(dex) 


M acc Reference 
(Moyr" 1 ) 


Orion Nebula Cluster -0.13 ± 0.02 Biazzo et al. (2011a) 

a Orionis -0.02 ± 0.09 Gonzalez-Hernandez et al. (2008) 


3 x 10" a Robbcrto ct al. (2004) 
3 x 10" 10 Ririiaco et al. (2011a) 


Taurus 0.00 ± 0.07 D'Orazi et al. (2011) 
Chamaeleon II -0.12 ± 0.14 This work 


3 x 10" 9 Giidel et al. (2007) 
7 x 10" 10 This work 



create planetesimals and protoplanets at the end of the 
CTTS phase. Protoplanets may clear the inner disk of 
gas and dust, causing the disk to lose its near- infrared 
color excess and at the same time opening a gap in the 
disk ([Lin fc Papaloizoulll993l) . thereby terminating accre- 
tion from the disk onto the star. 

With the aim of investigating possible relationships be- 
tween near-infrared colors and accretion properties, we used 
2MASS and SpitzerE colors (see Figs. [TO] and E]). We con- 
sidered these colors as disk tracers from the inner to the 
outer zone because they may estimate the magnitude of 
the near- and mid- infrared excesses above the photospheric 
level. In order to quantify the degree of correlation between 
these diagnostics, we calculated the Spearman's rank cor- 
relation coefficient, as we did in Section [4. 1[ and considered 
Class II stars with EWu a > 10 A. The correlation coeffi- 
cients, together with the probabilities, are listed in Tabled 
and show that the best agreements are obtained for the 
Ks— [8.0] and [3.6] — [4.5] colors versus (logM acc ) parame- 
ters. Also, Ks-[4.5], Ks-[5.8], [3.6]-[5.8], and [3.6] -[8.0] 
colors versus (logM acc ) show good agreement. This im- 
plies that objects with detectable accretion have optically 
thick inner disks. In particular, we can define the regions 
where M acc > 1.0 x 10" 10 M Q yr" 1 and K s-[8.0]>1.5 or 
[3.6]-[4.5]>0.2, or [3.6]-[5.8]>0.5, or [3.6]-[8.0]>0.9 as 
those where accreting objects with infrared excess are found 
in Cha II. The value M acc ~ lO _lo M yr~ 1 is a reasonable 
threshold for the t ransition from opticall y thick to optically 
thin (inner) d is k (iD'Alessio et al.|[2006t ). as also found by 
iRigliaco et all (|2011aD in the a Orionis SFR. The rough 
trend we tentatively observe among optically thin and op- 
tically thick disks, which needs to be confirmed on larger 
samples, might indicate a link between the mass accretion 
rate and the grain properties. This link, in turn controls the 
disk geometry, a connection that is worth exploring further. 

5. Conclusions 

In this paper, we determined radial velocities, lithium abun- 
dances, and accretion properties of 40 members of the 
Chamaeleon II star-forming region from FLAMES@VLT 
optical spectroscopy. Elemental abundances of Fe, Al, Si, 
Ca, Ti, and Ni for a suitable pre-main sequence star of the 
region were also measured. Our main results can be sum- 
marized as follows: 

1. The average radial velocity of the stars is consistent 
with that of the gas (~ 10 — 12 km s _1 ). The dispersion 
of the radial velocity distributions of the stars and gas 

5 We considered the IRAC@Spttzer fluxes at 3.6, 4.5, 5.8, and 
8.0 /an and the MlPS@Spitzer fluxes at 24 and 70 /im. 



Table 4. Spearman's correlation coefficients (p) and prob- 
abilities (p) for different color-M acc relations. 



Color 


P 


V 


Ks- 


[3.6] 


0.28 


24% 


Ks- 


[4.5] 


0.56 


5% 


Ks~ 


[5.8] 


0.52 


2% 


Ks- 


[8.0] 


0.68 


1% 


Ks- 


[24] 


0.29 


16% 


Ks- 


[70] 


0.30 


19% 


[3.6] 


-[4.5] 


0.80 


1% 


[3.6] 


-[5.8] 


0.61 


1% 


[3.6] 


-[8.0] 


0.52 


8% 


[5.8] 


-[8.0] 


-0.07 


80% 


[8.0] 


-[24] 


-0.09 


75% 


[24]- 


-[70] 


0.05 


85% 



(~ 1 — 2 km s 1 ) are also in agreem ent. Similar results 
were found by iDubath et al.l (1996) in Chamaeleon I 
(where (K a d)cha i ~ 15.0 ± 0.5 km s -1 ). 

2. A binary fraction of 18% is found in Chamaeleon II, 
which is similar to that in other T associations with 
comparable star density (see, e.g., the Taurus- Auriga 
association; iTorres et all 120021 ). 

3. The metallicity of the suitable member is slightly be low 
the solar value, as also found bv lSantos et al.l (|2008h for 
the Chamaeleon Complex. 

4. We find an average lithium abundance of the star- 
forming region of 2.5 ± 0.6 dex. 

5. Mass accretion rates derived through several secondary 
diagnostics (e.g., Ha, H/3, Hei A5876A, and Hei 
A6678A) are consistent with each other, justifying the 
use of all of them to compute an average mass accretion 
rate for each star. 

6. We provide a relationship between accretion luminosity, 
L acc , and the luminosity of the Hei A7065 A line. 

7. The relationship between mass accretion rate (M acc ) 
and stellar mass (M*) in Chamaeleon II is consistent 
with that found by previous studies in other T associa- 
tions. 

8. Although slightly lower, the average mass accretion rate 
in Chamaeleon I I fits well the relationship M acc versus 
Age reported bv ISicilia-Aguilar et al.l (|2010f ). 

9. We cannot exclude that significant variability on 
timescales longer than the time span of our observations 
and possibly due to episodes of variable mass accretion 
may produce the vertical scatter observed in the M acc 
versus M* plot. However, we suggest that such scatter 
of about two orders of magnitude in M acc at a given 
mass is also affected by other stellar properties. 
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Fig. 10. Infrared 2MASS-Spitzer colors versus mean mass accretion rates. Squares and dots correspond to Class III and 
Class II objects, respectively (see Table [7]). Symbol sizes represent stars with EW^a < 10 A (small squares and dots) 
and with EWn a > 10 A (big dots). Mean error bars are overplotted on the lower right corner of each panel. 



10. The color-M acc relationships suggest that the circum- 
stellar disks in Cha II become optically thin at ^ 
lO-^Moyr- 1 . 

11. The fraction of accretors in Chamaeleon II is ~ 50% 



in star-forming regions bv iFedele et all (12010ft . is con- 
sistent with the estimated age for the region (~ 3 Myr). 
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Table 5. Observing log, radial velocities, and lithium content. 



2MASS 
Name 




Other ID 


JD 
(+2 450 000) 


Instr. 


# 
obs. 


Kid 
(km s- 1 ) 


EW Li 

(A) 


(logn(Li) > 
(dex) 


Comment b,c 


12531722- 


-7707106 


IRAS 12496-7650/DK Cha 


3795.7806 


GIRAFFE 


2 









C/HV; SI, S2 






3797.7351 













C/HV; SI, S2 


12563364- 


-7645453 


Sz 46N 


3796.7383 


UVES 


3 


14.9±1.0 


0.61 


2.4 


Ol 








3797.8682 






15.9±0.5 


0.54 




Ol 








3797.8963 






15.8±2.0 


0.52 




Ol 


12565864- 


-7647067 


Sz 47 


3796.7382 


GIRAFFE 


3 









C/HV 








3797.8681 













C/HV 








3797.8963 


» 











C/HV 


12571172- 


-7640111 


IRAS 12535-7623/CHIIXR 2 


3796.7383 


UVES 


3 


10.2±1.3 


0.58 


2.9 










3797.8682 






9.9±1.5 


0.58 












3797.8963 


» 




9.9±1.8 


0.62 






12585614- 


-7630104 


WFI J12585611-7630105 


3796.7382 


GIRAFFE 


3 


11* 


0.53 


3.0 










3797.8681 






9* 


0.54 












3797.8963 






11* 


0.62 






12590984- 


-7651037 


C41 


3795.3476 


GIRAFFE 


1 




0.24 


1.8 


BC; SI, S2; Nl, N2 


12591013- 


-7712139 


Iso-Cha II 29 


4151.1828 


GIRAFFE 


2 




0.57 


2.6 


BC; SI, S2; Nl, N2 








4152.2866 








0.52 




BC 


13005346- 


-7709086 


Sz 48SW/CHIIXR 7 


3795.7394 


GIRAFFE 


3 


12* 


0.66 


2.6 








3796.8740 






12* 


0.63 












3797.6925 


1! 




11* 


0.63 






13005532- 


-7710222 


Sz 50/Iso-Cha II 52/CHIIXR 8 


3795.7394 


UVES 


3 


12.5±1.3 


0.52 


2.4 


Ol 






3796.8740 






13.9±1.2 


0.54 




Ol 








3797.6926 


»i 




12.6±1.1 


0.52 




Ol 


13005534- 


-7708296 


WFI J13005531-7708295 


3795.7394 


GIRAFFE 


1 


10* 


0.61 


2.6 




13005622- 


-7654021 


RX J1301.0-7654a 


3795.8653 


UVES 


1 


-4.2±2.4 


0.59 


3.7 


SB2 














32.6±2.7 








13005927- 


-7714027 


IRAS F12571-7657/Iso-Cha II 54 


3797.6925 


GIRAFFE 


1 


8* 


0.30 




LSN 


13015891- 


-7751218 


Sz 51/BC Cha 


3795.8943 


UVES 


2 


45.3±0.7 


0.42 


1.9 


SB1; Ol, 02; F 








3797.7786 






14.3±0.7 


0.30 




Ol, 02; F 


13021351- 


-7637577 


CM Cha/IRAS 12584-7621 


3795.8653 


GIRAFFE 


1 


13* 


0.43 


2.4 




13022287- 


-7734494 


C50 


3795.8936 


GIRAFFE 


2 


12* 


0.57 


2.5 










3797.7785 






13* 


0.41 






13030444- 


-7707027 


RX J1303. 1-7706 


3795.7394 


GIRAFFE 


3 


12.5±0.4 


0.61 


2.9 










3796.8740 






12.2±0.4 


0.62 












3797.6925 






12.3±0.4 


0.62 






13030905- 


-7755596 


C51 


3795.8936 


GIRAFFE 


2 


12* 


0.55 


3.1 










3797.7785 






11* 


0.56 






13042410- 


-7650012 


Hn 23 


3795.8653 


UVES 


1 


15.2±0.2 


0.53 


3.2 


Ol 


13042489- 


-7752303 


Sz 52 


3795.8936 


GIRAFFE 


2 


14* 


0.49 


1.9 


SI, S2; Nl, N2 








3797.7785 






13* 


0.48 




SI, S2; Nl, N2 


13045571- 


-7739495 


Hn 24 


3795.6959 


GIRAFFE 


3 


12±3 


0.60 


2.9 










3796.8263 






11±3 


0.61 












3797.6463 






11±3 


0.61 












3795.8943 


UVES 


2 


46.6±0.7 


0.61 


2.9 


SB1 








3797.7786 






14.1±0.8 


0.61 






13050855- 


-7733425 


Hn 25 


3795.6959 


GIRAFFE 


1 


10* 


0.47 


2.0 




13051269- 


-7730525 


Sz 53 


3795.6959 


GIRAFFE 


3 


8* 


0.40 


1.8 


SI, S2; Nl, N2 








3796.8263 






10* 


0.40 




SI, S2; Nl, N2 








3797.6464 






11* 


0.43 




SI, S2; Nl, N2 


13052072- 


-7739015 


Sz 54/BF Cha 


3795.8936 


GIRAFFE 


2 


14* 


0.46 


3.1 










3797.7785 






10* 


0.47 












3795.6960 


UVES 


3 


-1.0±1.4 


0.47 


3.1 


SB2; Ol 














30.9±1.6 














3796.8263 






-10.8±2.0 


0.47 




Ol 














34.0±1.6 














3797.6464 


» 




-6.3±2.3 


0.47 




Ol 














35.4±2.8 








13052169- 


-7738102 


SSTc2d J130521. 7-773810 


4151.2721 


GIRAFFE 


1 




0.33 




BC; C; SI, S2; Nl, N2 


13063053- 


7734001 


Sz 55 


3795.6959 


GIRAFFE 


3 


8* 


0.15 


0.6 


SI, S2; Nl, N2 








3796.8263 






13* 


0.26 




SI, S2; N2 








3797.6464 






13* 


0.27 




SI, S2; Nl, N2 


13063882- 


-7730352 


Sz 56 


3795.6960 


UVES° 


3 




0.56 


2.9 


very LSN; Ol 








3796.8263 


77 






0.50 




very LSN; Ol 








3797.6464 


77 






0.53 




very LSN; Ol 


13065656- 


-7723094 


Sz 57/C60 


3795.6960 


GIRAFFE 


3 


10* 


0.54 


3.0 








3796.8263 






10* 


0.57 












3797.6464 






10* 


0.53 






13065744- 


-7723415 


Sz 58/IRAS 13030-7707/C61 


3795.6960 


GIRAFFE 


3 


13±4 


0.49 


3.0 


SI, S2 






3796.8263 


77 




12±5 


0.49 




SI, S2 








3797.6464 


77 




12±4 


0.48 




SI, S2 


13070922- 


-7730304 


Sz 59/BK Cha 


3795.6959 


GIRAFFE 


3 


12±3 


0.39 


2.3 


S2 






3796.8263 


77 




14±3 


0.43 




S2 








3797.6464 


77 




12±2 


0.41 




S2 



Table 5. continued. 



2MASS 




Other ID 


JD 


Instr. 


# 


ya 
rad 


EW hi 


(logn(Li) > 


Comment b,c 


Name 






(+2 450 000) 




obs. 


(km s — 1 ) 


(A) 


(dex) 




13071806- 


-7740529 


C62 


3795.6959 
3797.6463 


GIRAFFE 


2 


13* 

8* 


0.65 
0.42 


3.0 




13072241- 


-7737225 


Sz 60W 


3795.6960 
3796.8263 
3797.6464 


UVES 


3 


13.9±0.8 
13.2±1.0 
12.5±0.9 


0.49 
0.51 
0.53 


2.2 




13074851- 


-7741214 


Hn 26 


3795.6959 
3796.8263 
3797.6463 


GIRAFFE 

" 


3 


11* 
10* 
10* 


0.59 
0.57 
0.60 


2.3 




13080628- 


-7755051 


Sz 61/BM Cha 


3796.6956 
3797.8201 


GIRAFFE 


2 


13±3 
11±3 


0.40 
0.39 


2.6 




13082714- 


-7743232 


C66 


3795.6959 
3796.6956 
3797.6463 
3797.8201 


GIRAFFE 

5J 
1) 
5J 


4 


9* 
8* 
11* 
12* 


0.55 
0.63 
0.45 
0.30 


2.5 


SB1? 


13090987- 


-7709437 


IRAS F13052-7653NW/CHIIXR 60 


4151.2219 


GIRAFFE 


1 




0.40 


1.9 


BC; S2 


13091071- 


-7709443 


IRAS F13052-7653N/CHIIXR 60 


3795.8228 


UVES 


1 


46.7±0.8** 


0.61 


2.5 


SB1? 


13095036- 


-7757240 


Sz 62 


3796.6956 
3797.8201 


GIRAFFE 


2 


10* 
10* 


0.44 
0.44 


1.8 




13100415- 


-7710447 


Sz 63 


3795.8227 


GIRAFFE 


1 


10* 


0.54 


2.2 




13125238- 


-7739182 


2MASS J13125238-7739182 


3796.6956 
3797.8201 


GIRAFFE 


2 


15* 
13* 


0.58 
0.60 


3.4 




13140369- 


-7753076 


Sz 64 


3796.6956 
3797.8201 


GIRAFFE 


2 


10* 
10* 


0.43 
0.43 


2.1 





Notes: 

* Due to low S/N ratio, few lines, late spectral type, and/or short wavelength coverage, the radial velocity error may be up to 
60%. 

** This can be classified as a bona-fide PMS star based on other criteria (see Alcala ct al. 2008; Spezzi et al.ll2008D ; its high 
radial velocity suggests it may be a SB1. 

a For spectroscopic binaries, the different radial velocity measurements of the corresponding CCF peaks are given. 

b SB1: single-lined spectroscopic binary; SB2: double-lined spectroscopic binary; SB1?: suspected SB1; BC: bad CCF; C: 

continuum-type spectrum; HV: heavily veiled star; LSN: low S/N spectrum. 

c 01=[Oi] A6300.8A, 02=[Oi] A6363.8A, Sl=[Sli] A6715.8A, S2=[Sn] A6729.8A, Nl=[Nn] A6548.4A, N2=[Nn] A6583.4A, 
F= [ Fen] multiple ts. 

° In lSpezzi etaLI (|2008T ). it was mistakenly reported as GIRAFFE observations. 



Table 6. Star name (columnl); equivalent width, observed flux, and mass accretion rate of the Ha, H/3, He I A5876 A, and A6678 
A lines (columns 2-13); He I A 7065 A line equivalent width and observed line luminosity (columns 14-15); average M acc as derived 
from the Ha, H/3, He I A5876 A, and A6678 A lines (column 16). 



ID narnc 


EW Ha 

( A 1 
(A) 


log F H <* log M"° 
(ergs -1 cm 2 ) (M yr" 1 ) 


EW H p 

( A "i 
(A) 


log F Hp 

^crg s cm ) 


loe M"' EWh,.™-tr 

& acc L ''^ilcoo7o 

(M & yr ) (A) 


log F Hc ^ 7H 

[erg s cm ) 


(M Q yr~ ) 


/' l V U" & ^ 7C 

y ' .ricOOYO 


log F n 
(ergs -1 


eS678 
cm ) 


log M£j ebYS 
(Mq yr -1 ) 


Hi VV Hc7065 


logL Hc70 ^ 


\ 1U 6 Jvl acc / 
{Mq yr ) 


IRAS 12496 — 7650/DK Cha 


78.5 


9.1 


-6.3 












0.65 




7.0 


— 6.1 






— 6.1 


91.5 


9.2 


-6.2 












1.07 




7.2 


— 5.8 








Sz 46N 


10.5* 


6.7 


-9.7 


4.6 


6.1 


— 9.8 0.50 


5.2 


— 9.6 


0.06 




4.2 


— 10.5 






— 9.5 




15.6* 


6.8 


-9.5 


8.5 


6.4 


— 9.5 0.65 


5.3 


— 9.4 


0.07 




4.3 


— 10.4 










17.4* 


6.9 


-9.4 


5.9* 


6.2 


— 9.7 0.69 


5.3 


— 9.4 
















Sz 47 


32.5 
































32.5 
































32.3 






























T"r> AO ioror ^7/^00 //^ttttv t> 

IKAb 12535 — 7623/CIIIIJyK. 2 


2.2* 


6.0 


-9.6 


0.8 


5.6 


— 9.6 


















— 9.8 




1.5 


5.9 


-9.8 


0.6* 


5.5 


— 9.7 






















l.o 


5.8 


-9.8 


0.3* 


5.2 


— 10.1 




















WFI J12585611 — 7630105 


29.2 


6.3 


-9.7 












0.35 




4.2 


— 10.1 






— 9.9 




31.2 


6.3 


-9.7 












0.55 




4.4 


— 9.8 










31.5 


6.3 


-9.7 












0.52 




4.4 


— 9.9 








C41 


46.7 


6.8 


-10.8 












0.95 




4.9 


— 10.9 


0.18 


— 7.3 


— 11.0 


Iso-Ona 11 29 


0.9 


5.6 


-10.6 
























— 10.6 




1.0 


5.7 


-10.5 


























bz 48bW/CHIIXK 7 


9.9* 


6.6 


-9.9 
























— 9.9 




14.7* 


6.8 


-9.7 




























8.1* 


6.6 


-10.0 


























O m /T TT [~ O /riTTTTVD o 

bz 50/lso-Cna 11 52/CHIIJyK 8 


33.8* 


6.8 


-8.1 


9.7 


5.9 


— 8.6 0.43 


4.6 


— 8.7 


0.09 




4.1 


— 9.1 






— 8.4 




35.4* 


6.8 


-8.1 


11.8* 


6.0 


— 8.5 0.46 


4.7 


— 8.7 


0.08 




4.0 


— 9.2 










30.4* 


6.7 


-8.1 


8.5 


5.9 


— 8.6 0.32 


4.5 


— 8.9 


0.07 




4.0 


— 9.2 








WFI J13005531 — 7708295 


2.1 


6.0 


-10.4 
























— 10.4 


Ha J1301.0 — 7654a 


3.1 


6.8 


-8.6 


1.4 


6.5 


— 8.4 0.16 


5.6 


— 7.9 














— 8.3 


IHAb 1 12571 — 7657/Iso-Cna II 54 


16.0 






























Sz 51 


132.5* 


8.1 


-8.1 


39.0 


7.5 


— 8.3 2.92 


6.4 


— 8.0 


1.11 




6.0 


— 8.3 






— 8.3 




110.0* 


8.0 


-8.2 


26.0* 


7.3 


— 8.5 2.29 


6.3 


— 8.1 


0.91 




5.9 


— 8.4 








CM Cna/lKAb 12584 — 7621 


27.6 


7.4 


-8.6 












0.72 




5.8 


— 8.2 


0.50 


— 4.6 


— 8.4 


C50 


37.5 


6.6 


-9.8 












1.34 




5.1 


— 9.5 


0.39 


— 6.2 


— 9.9 




38.3 


6.7 


-9.8 












0.78 




4.8 


— 9.9 


0.48 


— 6.1 




Ha J1303. 1 — 7706 


2.6 


6.1 


-9.5 
























— 9.5 




2.7 


6.1 


-9.5 




























2.6 


6.1 


-9.5 


























C51 


9.5 


6.0 


-9.6 
























— 9.7 




6.3 


5.8 


-9.8 


























Hn 23 


11.9 


7.4 


-8.9 


1.5 


6.5 


— 9.5 0.26 


5.7 


— 8.8 


0.08 




5.2 


— 9.2 






— 9.0 


Sz 52 


49.4 


7.1 


-9.2 












0.70 




5.1 


— 9.4 


0.45 


— 5.6 


— 9.4 




49.9 


7.1 


-9.2 












0.62 




5.1 


— 9.4 


0.43 


— 5.6 




Hn 24 


1 Q« 

1.3 


5.8 


-10.0 
























— 9.8 




1.2 


5.8 


-10.0 




























1.9 


6.0 


-9.7 




























1.1 


5.7 


-10.0 


0.9 


5.7 


— 9.6 0.12 


4.6 


— 9.5 


















1.4 


5.8 


-9.9 


0.8 


5.6 


— 9.7 




















Hn 25 


22.5 


6.6 


-9.4 












0.67 




4.9 


— 9.2 


0.54 


— 5.4 


— 9.3 


Sz 53 


145.0* 


7.8 


-8.2 












1.05 




5.5 


— 8.7 


1.37 


— 4.7 


— 8.6 




46.9* 


7.3 


-8.8 












0.70 




5.3 


— 8.9 


0.64 


— 51 






98.3* 


7.6 


-8.4 












0.91 




5.4 


-8.8 


0.79 


-5.0 




Sz 54 


23.8 


7.7 


-7.8 
























-8.2 




22.7 


7.7 


-7.8 




























21.2 


7.7 


-7.9 


1.9 


6.6 


-8.6 






















22.4 


7.7 


-7.8 


2.2 


6.7 


-8.5 






















22.6 


7.7 


-7.8 


1.2 


6.4 


-8.8 




















SSTc2d J130521. 7-773810 


26.1 
















0.21 














Sz 55 


101.1* 


7.6 


-9.1 












2.04 




5.7 


-9.1 


1.39 


-5.3 


-9.1 




127.2* 


7.7 


-9.0 












1.93 




5.7 


-9.1 


1.41 


-5.2 






130.8* 


7.7 


-9.0 












1.97 




5.7 


-9.1 


1.14 


-5.3 





Table 6. continued. 



ID name 


EW Ha 


logF Ha 




EWnp 


log F Hf * 


log AfW 




logi rHe5S7<3 


logM a H c f™ 


Whc6678 


lo£ 


. pHeGGTH 


log M a H c f 6 ™ 




logL"^ 


(logM acc ) 




(A) 


(ergs -1 cm 2 ) 


(Mgyr- 1 ) 


(A) 


(ergs -1 cm 2 ) 


(M yr- 1 ) 


(A) 


(erg s - 1 cm 2 ) 


(Mgyr- 1 ) 


(A) 


(erg 


s-'em 2 ) 


(Mgyr -1 ) 


(A) 


(£©) 


(Mgyr -1 ) 


Sz 56 


11.1* 


6.1 


-9.7 


3.5* 


5.1 


-10.3 


1.90 


5.0 


-9.0 


0.60 




4.7 


-9.1 






-9.4 




14.3* 


6.2 


-9.5 


12.8* 


5.7 


-9.7 


3.70 


5.3 


— 8.6 


0.78 




4.8 


-8.9 










18.8* 


6.3 


-9.4 


4.8* 


5.2 


-10.2 


2.97 


5.2 


-8.8 


0.88 




4.8 


-8.9 








Sz 57/C60 


23.0* 


6.4 


-8.7 














0.50 




4.6 


-8.5 


0.12 


-5.9 


-8.7 




16.5* 


6.3 


-8.8 














0.39 




4.5 


-8.7 


0.11 


-5.9 






16.2* 


6.3 


-8.9 














0.41 




4.5 


-8.7 


0.10 


-6.0 




Sz 58/IRAS 13030-7707/C61 


18.2* 


7.6 


-8.8 






















0.15 


-5.0 


-8.9 




14.8* 


7.5 


-8.9 














0.02 




4.6 


-10.2 


0.21 


-4.9 






10.5* 


7.3 


-9.1 














0.08 




5.2 


-9.4 


0.20 


-4.9 




Sz 59 


52.8* 


7.7 


-8.1 














0.44 




5.6 


-8.3 


0.48 


-4.5 


-8.4 




35.3* 


7.5 


-8.3 














0.30 




5.4 


-8.5 


0.30 


-4.7 






39.1* 


7.6 


-8.3 














0.25 




5.3 


-8.6 


0.30 


-4.7 




C62 


37.2 


6.6 


-9.7 














1.40 




5.0 


-9.4 


0.45 


-6.1 


-9.6 




39.6 


6.6 


-9.7 














1.85 




5.2 


-9.2 


0.60 


-6.0 




Sz 60W 


25.6* 


7.1 


-9.3 


7.5* 


6.4 


-9.6 


1.22 


5.5 


-9.1 


0.38 




5.0 


-9.5 






-9.3 




37.6* 


7.2 


-9.1 


11.8* 


6.5 


-9.4 


0.93 


5.4 


-9.3 


0.31 




4.9 


-9.6 










36.8* 


7.2 


-9.1 


12.4* 


6.6 


-9.4 


1.06 


5.5 


-9.2 


0.39 




5.0 


-9.5 








Hn 26 


6.6* 


6.3 


-10.1 














0.08 




4.3 


-10.5 


0.09 


-6.1 


-10.1 




8.1* 


6.5 


-10.0 














0.18 




4.6 


— 10.0 


0.08 


— 6.1 






13.4* 


6.7 


-9.7 














0.11 




4.4 


-10.2 


0.11 


-6.0 




Sz 61/BM Cha 


81.1* 


8.2 


-7.6 














0.68 




6.1 


-7.7 


0.80 


-4.0 


-7.7 


64.7* 


8.1 


-7.8 














0.60 




6.1 


-7.7 


0.53 


-4.2 




C66 


20.8 


6.5 


-10.2 














0.67 




4.9 


-10.0 


0.50 


-6.2 


-9.8 




29.1 


6.7 


-10.0 














1.87 




5.3 


-9.4 


0.40 


-6.3 






35.9 


6.8 


-9.9 














2.03 




5.4 


-9.3 


0.99 


-5.9 






25.7 


6.6 


-10.0 














2.10 




5.4 


-9.3 


0.90 


-5.9 




IRAS F13052-7653NW/CHIIXR 60 


11.2 


6.8 


-10.0 














0.73 




5.4 


-9.4 


0.50 


-5.3 


-9.7 


IRAS F13052-7653N/CHIIXR 60 


2.2 


5.9 


-10.5 


0.5 


5.0 


-11.0 




















— 10.8 


Sz 62 


115.8 
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-8.4 
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0.98 




5.0 
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0.88 


-5.2 




Sz 63 


93.6 


7.2 


-8.9 














0.41 




4.6 


-9.8 


0.35 


-5.8 


-9.4 


2MASS J13125238-7739182 


5.8 


5.8 


-10.3 














0.11 




3.9 


-10.5 






-10.5 




7.0 


5.9 


-10.2 














0.11 




3.9 


-10.5 








Sz 64 


132.1 


7.1 


-9.1 














2.82 




5.3 


-9.0 


0.79 


-5.9 


-9.1 




121.4 


7.0 


-9.1 














3.00 




5.3 


-9.0 


1.22 


-5.7 





Notes: 

— * Variable line profile. 
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Notes: 

a Spectral types, effective te mperatures, l u minos ities, and radii are from Spe zzi et alJ l|2008l '). Masses and ages are the average 
of three values reported in ISpezzi et al.l (|200Sf ) and derived using three different sets of evolutionary models, with mean 
standard deviations of ~ O.15M0 and ~ 1.8 Myr, respectively. Mean surfa ce gravities were co mputed through the relation: 
(loggr) = 4.4377 + (log M+) + 4 log T eff /5777 - logL*. R in was adopted from lAlcala et al.1 (|2008l ). 
b Lada Class as derived in lAlcala et al.l l |2008f l. 

* In lAlcala et"a l. (2008), it appears as Class III, but it is a Class II object. 

* Rin = 5R+ (see Section EH). 
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Appendix A: Examples of Ha/H/3 line profiles 

Here, we display the Ha (Figs. IA.11 IA.2|) and H/3 (Figs. IA.3|) profiles of the targets showing changes in their line shape 
and/or intensity. Below, we briefly comment on each object. 

Sz 48SW. During the first observing night, the spectrum showed a symmetric, relatively narrow Ha emission profile with 
a peak close to the line center (type I profile, following the classification of iReipurth et ai1ll996f) . During the second 
night, the profile appeared slightly asymmetric with lower emission on the blue than on the red side (type HIB), and 
then turned back to the type I profile on the third night. 

Hn 24- The Ha line profile changes from type IVR (an inverse P-Cygni-like profile) during the first two nights to type 
IIIR, where less emission in the red than in the blue is seen. 

Sz 53. The Ha line shows always a type IIB profile, with a central reversal at the line center and the blue-wing peak 
intensity lower than the red one. The line strength changes, being maximum during the first observation. The secondary 
peak always exceeds half the strength of the primary peak. This type of peak is probably due to the interplay of variable 
accretion and mass loss. 

Sz 55. This star shows an Ha line profile changing from IIB to IIIB (where the secondary peak is slightly less intense 
than half the strength of the primary peak) to IIB. 

Sz 57. It always shows a narrow and symmetric Ha emission type I profile with a peak close to the line center. 

Sz 58. The Ha profile changes from IIIR (first night) to IIR (second and third night). 

Sz 59. IIR/IIR/IIIR profiles are observed for this star during the three observations. 

Hn 26. This star shows IIR/IIB/IIB Ha line profiles. 

Sz 61. It shows an Ha profile slightly changing from IIB type to IIIB. 

Sz 46N. This star shows Ha and H/3 line profiles that are always symmetric (type I). 

IRAS 12535-7623. The Ha and H/3 emission line profiles are always consistent with type I. 

Sz 50. Its Ha and H/3 emission line profiles are always narrow and almost symmetric (type I). 

Sz 51. Ha and H/3 are always of type I. 

Sz 56. The Ha profile changes from IIR (during the first and second night) to IIIR (third night), while the H/3 profile is 
always of type I. 

Sz 60W. The Ha line has a complex profile, starting from IIB type during the first night; during the second night, it 
shows broad wings, red-shifted and blue-shifted absorption, and a narrow emission in the center (like a YY Orion Ha 
profile, normally associated with high infall and outflow rates, consistent with the high value of M acc ; IWalkerl Il972h : 
during the third night, it turned back to a IIB type profile. The H/3 line changes from YY Orionis H/3-likc profile to type 
IIIR to type I. 



Appendix B: Stellar variability 
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Fig. A.l. Examples of Ha emission line profile variations of nine stars. The fluxes are normalized to the continuum. The 
data refer to the FLAMES/GIRAFFE configuration. The solid/dashed/dash-dotted line represents the first/second/third 
observation, respectively. 




Fig. A.2. Same as in Fig. [O but for the UVES data. 
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Fig. A. 3. Same as in Fig. [O but for the UVES data and the R/3 profiles. 
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Fig. B.l. Mass accretion rates measured from the Ha (left top panel), H/3 (right top panel), Hei A5876A (left bottom 
panel), and He I A6678 A (right bottom panel) lines versus mass. Vertical bars connect the minimum and maximum values 
of Mace obtained per each star at different epochs. Stars are divided into high Ha emitters (EWu a > 10 A; grey bars) and 
low Ha emitters (EW-aa < 10 A; orange bars), i.e., stars where the line emission cannot be unambiguously attributed to 
accretion activity, as it is most likely caused by chromospheric activity. The dashed line shows the relation M acc oc M^- 87 
obtained by iHerczeg fc Hillenbrand! (|2008l ) for Taurus members. The highest M acc values in the first and fourth panels 
refer to the early-type star DK Cha. The lower number of points in the H/3 and He I plots is due to two reasons: absence 
of the line in the given star or line out of the wavelength range (see Table [5]). 



